RESEARCH PAPER 



Plant Signaling & Behavior 8:8, e24937; August 201 3; © 201 3 Landes Bioscience 

Herbivory of maize by southern corn rootworm 
induces expression of the major intrinsic protein 
ZmNIPI;! and leads to the discovery of a novel 

aquaporin ZmPIP2;8 

Susan D. Lawrence/* Nicole G. Novak/ Hao Xu^ and Janice E.K. Cool<e^ 

'USDA-ARS; Invasive Insect Biocontrol and Behavior Lab; BARC-West; Beltsville, MD USA; ^Department of Biological Sciences; University of Alberta; Edmonton, AB Canada 
Keywords: aquaporin, corn rootworm, Diabrotica, herbivory, maize, root 

Abbreviations: SCR, southern corn rootworm; MIP, Major intrinsic protein; NPA, Asparagine-Prohne-Alanine; ar/R, aromatic/ 
arginine; NIP, Nod-26-hke intrinsic protein; PIP, plasma membrane intrinsic protein; TIP, tonoplast intrinsic protein; TM, 
transmembrane; RIP2, ribosome-inactivating protein; SDP, specificity determining residues 



Aquaporins channel water and other neutral molecules through cell membranes. Aquaporin gene expression is subject 
to transcriptional control and can be modulated by factors affecting water balance such as salt, abscisic acid and drought. 
During infestation of maize by southern corn rootworm (SCR), an insect that chews into and significantly damages maize 
roots, three maizeaquaporins were differentially expressed upon prolonged infestation. Using a brief infestation of maize 
roots ZmA//P7;7 transcript abundance again increased under infestation while expression of a new aquaporin, ZmP/P2;S and 
ZmTIP2;2 expression did not change. Since ZmPIP2;8 has not been described previously, the deduced protein sequence 
was analyzed in silico and found to contain the hallmarks of plant aquaporins, with a predicted protein structure similar 
to other functionally characterized PIP2s. NIPs characterized to date have been implicated in facilitating the movement 
of a variety of small molecules, while TIPs and PIPs often have the capacity to facilitate trans-membrane movement of 
water. Functional assays (using heterologous expression in Xenopus laevis oocytes) of ZmTIP2;2 and ZmPIP2;8 confirmed 
that these aquaporins demonstrate water channel capacity. 



There are at least one million species of phytophagous insects 
that rely on plants for their survival, so it is not surprising that 
plants have evolved well-orchestrated responses to herbivory. 
Because they are sessile organisms, plants must utilize rapid and 
efficient signal transduction strategies to cope with predation 
and ensure survival. Many microarray studies have illustrated the 
extensive transcriptional reprogramming a plant may experience 
while under insect attack.''^ Whole genome shifts in transcript 
production in response to herbivory are thought to represent 
reprioritization to balance basic physiological processes such as 
growth and reproduction with general stress related responses.''''^ 
Root feeders present particular challenges to the plant, as 
they destroy tissues involved in nutrient acquisition, storage and 
water balance; processes that are all essential for plant survival. 
Therefore, many plant responses to this underground herbivory 
are not only defense strategies (i.e., chemical deterrents, vola- 
tile emission etc.), but also those that affect primary metabo- 
lism, resource allocation, osmotic regulation and growth.^ It is 
in osmotic regulation and growth where the ubiquitous integral 



membrane proteins known as Major Intrinsic Proteins (MIPs) 
figure. Aquaporins, the water transporting MIPs that facilitate 
movement of water across membranes, play key roles in osmotic 
regulation and growth-driving cell expansion, for review see ref- 
erences 8-9. 

Previously, a whole genome microarray was utilized to elu- 
cidate the effects of severe southern corn rootworm Diabrotica 
undecimpunctata howardi (Barber) (SCR) infestation on the 
maize transcriptome.^ In this study, differentially expressed 
MIPs were annotated using the Maize Genome Database (www. 
maizegdb.org). It was found that transcript abundance cor- 
responding to ZmNIPl;! was upregulated under conditions of 
infestation while that of ZmTIP2;2 and a new, previously unre- 
ported ZmPIP2 were downregulated. 

MIPs are found in virtually every species of mammals, plants 
and microorganisms.'* ''' While the water-transporting aquaporins 
play a crucial role in water relations, the transport of non-water 
substrates such as ammonia, hydrogen peroxide and glycerol by 
some members of the MIP family suggests a larger role in plant 
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physiology.' An extensive overview of 470,000 EST sequences 
representing 215 maize cDNA libraries identified 33 unique 
aquaporins." Phylogenetic analysis of the maize aquaporins 
identified four sub-families; plasma membrane intrinsic proteins 
(PIPs), tonoplast membrane intrinsic proteins (TIPs), small basic 
intrinsic proteins and nodulin26-like intrinsic proteins (NIPs). 
The term NIP was coined from the first identified NIP, soybean 
Nod-26, which is involved in a symbiotic interaction in the peri- 
bacteroid membrane of root nodules.' The X intrinsic proteins, 
another subfamily are found in several plant taxa, but have not 
so far been identified in some major plant groups such as the 
Brassicaceae or monocots.'^'''' 

Structurally, all MIPs have six trans-membrane helices sepa- 
rated by five loops; three extracellular (Loop 1, 3 and 5) and two 
intracellular (Loop 2 and 4). The helices orient longitudinally as 
a tetramer in the membrane to form a primary pore skeleton that 
substrates travel through. Inside this primary scaffold, two key 
areas or "constrictions" contribute to the pore architecture and 
substrate transport specificity. The asparagine-proline-alanine 
(NPA) residues in Loop 2 and 5 are hydrophilic and embed into 
the membrane to form one constriction. Selectivity for water is 
important in this area. The NPA constriction also produces a 
strong electrostatic barrier, effectively inhibiting the passage of 
protons and maintaining the important electrochemical gradi- 
ent across the membrane.**'" The second, narrower constriction is 
known as the aromatic/arginine (ar/R) constriction or filter and 
is formed inside the molecule from only four residues; two from 
transmembrane (TM) helix 2 and two from Loop 5. Here sub- 
strate passage is limited by size and hydrophilic or hydrophobic 
character of the residues in the constriction. It is also a secondary 
site of proton exclusion.' Additional conserved amino acids also 
influence whether glycerol or water are transported. Froger et al."" 
identified five additional conserved amino acids P1-P5 by exam- 
ining protein sequences of 40 aquaporins, which predict whether 
the MIP transports water or glycerol. PI is found in loop C, P2— 3 
within loop E and P4— 5 within TM6. The P1-P5 consensus sig- 
nature of glycerol channels is (Y/F)D(K/R)P(L/M/I/A/V) while 
the water channels are (Q/M/E/T/A)S(S/A)(F/Y)W. 

While transcriptional regulation of aquaporins in plants by 
abiotic stress has been well established,' less is known about tran- 
scriptional regulation of aquaporins by biotic stresses, particu- 
larly insect herbivory. Here, we assess transcriptional response of 
three aquaporins by SCR herbivory, one of which has not been 
previously described. In silico analyses were used to identify con- 
served aquaporin motifs in the sequence of the new PIP2, des- 
ignated ZmPIP2;8, and Xenopus oocyte swelling assays used to 
assess the water channel function of ZmPIP2;8 and ZmTIP2;2. 

SCR Infestation Response of Maize Aquaporin 
Transcripts 

Previous microarray results demonstrated that infestation for 
four weeks induced ZmNIPl;! and repressed ZmTIP2;2 and 
ZmPIP27 Infestations were performed on 7-d-old maize plants 
with six second to third instar larvae for 2 to 4 d. Using qRT- 
PCR, the three aquaporin genes were assayed (Fig. 1) as was a 



ribosome inactivating protein {ZmRIP2) . Although ZmRIP2 is 
not a MIP it has been previously characterized as an insect defense 
gene in maize and therefore used to confirm a defense response to 
this brief infestation.'^ ZmRIP2 and ZmNIPl;! transcripts were 
induced by infestation. However the ZmTIP2;2 and ZmPIP2;8 
were not differentially expressed by this brief infestation. 

ZmPIP2;8 Amino Acid Sequence Structure 

The ZmPIP2 sequence was predicted to have six a helices and 
5 loops with loop B and D inside and loops A, C and E out- 
side the cell. Both trans-membrane domain-predicting programs 
generate similar proteins, confirming the structure of a typical 
PIP. However each program places the helices at slightly differ- 
ent locations (Table 1). An amino acid sequence alignment of 
all 8 ZmPIP2s is shown (Fig. 2) and highlights the conserved 
NPA domains of loop B and E, which are thought to form the 
pore" along with an invariant H residue in loop D, the Ar/R 
filter is also shown. The TM domains are also marked (Fig. 2), 
using information from Hove and Bhave"* and the results from 
the TM domain programs as a guide. The similarity of this novel 
PIP to the other ZmPIP2s varies from 83—87%. Comparison to 
ZmPIPl proteins show at best 80% similarity to ZmPIPl;3/l;4. 
Since seven ZmPIP2 proteins have been described in maize," this 
283 amino acid protein was named ZmPIP2;8. A phylogenetic 
tree comparing ZmPIP2;8 to other PIP2 proteins from maize 
was constructed (Fig. 3). While ZmPIP2;8 is most closely associ- 
ated with ZmPIP2;2 it is clearly a unique ZmPIP2. 

Functional Assay of ZmTIP2;2 and ZmPIP2;8 

The ZmTIP2;2 and ZmPIP2;8 clones were assayed for their pro- 
tein's capacity to channel water using cRNA injection oi Xenopus 
oocytes as described in the "Materials and Methods." Neither of 
these proteins has been tested previously for this function. The 
results were compared with the strong water permeability of 
oocytes injected with cRNA of L. bicolor aquaporin JQ585595 
as the positive control (Fig. 4). The values for oocytes micro- 
injected with cRNAs of ZmPIP2;8 and ZmTIP2;2 were respec- 
tively, 11- and 6.1-fold higher than those of oocytes injected 
with water, showing ZmPIP2;8 and ZmTIP2;2 were functional 
aquaporins with significant water transport capacity in this het- 
erogonous expression system (Fig. 4). The Pj. value of the protein 
corresponding to ZmPIP2;8 was significantly higher than that 
corresponding to ZmTIP2;2. 

In experiments for the microarray, roots were exposed to SCR 
larvae for a longer period of time. Wilting and browning of those 
plant shoots may indicate a downregulation of water transport- 
ers such as ZmPIP2;8 and ZmTIP2;2, which can significantly 
change hydraulic conductivity.' 

Brief infestation does not have a similar effect on these trans- 
porters nor does it result in any obvious shoot stress. Further 
studies are required to correlate the role of infestation on the 
water transporters ZmTIP2;2 and ZmPIP2;8. The induction 
of ZmNIPl;! under both a brief and prolonged infestation may 
suggest that it plays a different role during herbivory other than 
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Figure 1. SCR infestation of maize roots significantly induced ZmNIPI;! and ZmRIP2 (AAF60304) but did not affect ZmPIP2;8 and ZmTIP2;2. Plants were 
grown, harvested and RNA was isolated from maize roots. cDNA was synthesized and qRT-PCR performed as described in Lawrence et a\7 Relative 
Quantity (RQ) was calculated where control uninfested RNA at time 0 was designated as the calibrator and scaled to a value of I.The data are ex- 
pressed as \og^^; the log^^^ transformed value of the control is zero. Error bars represent standard deviation. Three biological replicates are averaged for 
each time point. Four plants were pooled for each biological replicate. 



responding to changes in osmotic balance. Perhaps this MIP 
is transporting an alternative substrate involved in the plant's 
defense response. 

The conserved amino acids of the MIPs affected by SCR 
infestation were examined. Considering that the ZmTIP2;2 and 
the ZmPIP2;8 transport water, do these proteins contain the 
conserved amino acids from the Ar/R filter or P1-P5 residues 
with the properties of a water channeling aquaporin? A compari- 
son of conserved amino acids of the MIPs to AQPl the proto- 
typical water channel protein and GlpF the prototypical glycerol 
transporter (Table 2) shows that in the Ar/R filter and the P1-P5 
residues of ZmPIP2;8 and ZmTIP2;2 share six and seven of the 
nine residues respectively with the water channel protein AQPl. 
Although the Ar/R filter contains two nonpolar neutral residues 
at the second and third position similar to GlpF, all P1-P5 resi- 
dues of ZmTIP2;2 are common to water channel proteins such 
as AQPl."^ The ZmPIP2;8 differs from AQPl only at the third 
position of the Ar/R filter, which is also a T for all the other 
ZmPIP2 proteins and the PI residue has an H rather than T for 
AQPl or an invariant Qfor the other ZmPIP2s. This last differ- 
ence is a change from what is a polar neutral amino acid for all 
the ZmPIP2s to a basic residue for ZmPIP2;8. These differences 
from AQPl however still allowed functional water transport in 
Xenopus oocytes (Fig. 4). 



Table 1. Prediction of secondary structure of ZmPIP2;8 according to 
TMHIVIIVI version 2.0 and HIVIMTOP version 2.0 
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ZmPIP2;8 is 283 amino acids and the numbers represent the position of 
the protein sequence within each domain as predicted by the different 
programs. Both programs predict six trans-membrane helices and that 
the N-terminal sequence and the C-terminal sequence reside inside 
(intracytoplasmic) the cytoplasm along with loop B and D, while loop A, 
C and E are outside (extracytoplasmic) the membrane. 
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ZmPIP2;1 
ZmPIP2;2 
ZmPIP2;3 
ZmPIP2;4 
ZmPIP2;5 
ZmPIP2;6 
ZmPIP2;7 
ZmPIP2;8 



ZmPIP2;1 
ZmPIP2;2 
ZmPIP2;3 
ZmPIP2;4 
ZmPIP2;5 
ZmPIP2;6 
ZmPIP2;7 
ZmPIP2;8 



ZmPIP2;1 
ZmPIP2;2 
ZmPIP2;3 
ZmPIP2;4 
ZmPIP2;5 
ZmPIP2;6 
ZmPIP2;7 
ZmPIP2;8 



ZmPIP2;1 
ZmPIP2;2 
ZmPIP2;3 
ZmPIP2;4 
ZmPIP2;5 
ZmPIP2;6 
ZmPIP2;7 
ZmPIP2;8 



ZmPIP2;1 
ZmPIP2;2 
ZmPIP2;3 
ZmPIP2;4 
ZmPIP2;5 
ZmPIP2;6 
ZmPIP2;7 
ZmPIP2;8 

ZmPIP2;1 
ZmPIP2;2 
ZmPIP2;3 
ZmPIP2;4 
ZmPIP2;5 
ZmPIP2;6 
ZmPIP2;7 
ZmPIP2;8 



MGKDDVI ESGAGGGEFAAKDYT 
MGKDDVVQSGAGGGEFAAKDYT 
MAKQDI EASGPEAGEFSAKDYT 
MAK- Dl EASGPEAGEFSAKDYT 

MAK- Dl EAAAAHE GKDYS 

MGK- EVDVSTL EAGGVRDRDYA 
MAK- DVEQVTEQG- EYSAKDYH 
MDTCHQSI ETADG- - - - KKDYS 

TMJ 

TLLFLYI TVATVI GYKHQT 
TLLFLYVTVATVI GYKHQT 
TLLFLYI TVATVI GYKHQT 
TLLFLYI TVATVI GYKHQT 
TLLFLYI TVATVI GYKHQT 
TLLFLYI TVATVI GYKHQT 
TLLFLYI TVLTI I GYKRQS 
TLLFVYVTLATVI GHKREA 



DPPPAPL 
DPPPAPL 
DPPPAPL 
DPPPAPL 
DPPPAPL 
DPPPAPL 
DPPPAPL 
DPVPAPF 



I DAAELGSWSL 
VDAAELGSWSL 
I DADELTKWSL 
I DAEELTQWSL 
VDAEELTKWSL 
I Dl DELGKWSL 
I DPDELTKWSL 
VNAGELGKWSL 



YRAVI 
YRAVI 
YRAVI 
YRAVI 
YRAVI 
YRAVI 
YRAAI 
YRAVI 

TM2 



AEFI A 50 

AEFI A 50 

AEFI A 50 

AEFI A 49 

AEFVA 45 

AEFVA 49 

AEFI A 48 

AEFVA 46 



DASAS- - GADAACGGVGVLGI 
DASASGAGADAACGGVGVLGI 
DAAAS- - GPDAAGGGVGI LGI 
DASAS- - GPDAAGGGVGI LGI 
DAAAS- - GPDAACGGVGVLGI 
DASAS- - GPDAACSGVGI LGI 
DTKI P- - G- NTECDGVGI LGI 



ESQP CGSVGVLGI 



AWAFGGMI 
AWAFGGMI 
AWAFGGMI 
AWAFGGMI 
AWAFGGMI 
AWAFGGMI 
AWAFGGMI 
AWSFGGMI 



FV 

FV 

Fl 

Fl 

Fl 

Fl 

Fl 

FV 



TM2 



NPA LoopB 



LVYCT 
LVYCT 
LVYCT 
LVYCT 
LVYCT 
LVYCT 
LVYCT 
LVYCI 



AG 
AG 
AG 



SGGHI NPAVT FGI Fl ARKVRI 
FLARKVSL 
FL ARKVSL 
FLARKVSL 
FLARKVSL 
FLARKVSL 
FLGRKVSL 
LLARKLSL 



SGGHI NPAVT FGI 
RGGHI NPAVT FGI 
AGI SGGHI NPAVTFGI 
AGVSGGHI NPAVTFGL 
AGI SGGHI NPAVT FGL 
AGI SGGHI NPAVT FGI 
AGI SGGHI NPAVT FGL 



TM3 

vrALLVI VASCL5AI CfiVSLV 

VRALL YMVAQCLGAVCGVGL V 
VRALLYI I AQCLGAI CGVGLV 
VRALLYI I AQCLGAI CGVGLV 
VRALLYI VAQCLGAI CGVGLV 
VRALLYMAAQSLGAI CGVALV 
VRALLYMI AQCAGAI CGAGLA 
VRAAL YVVAQCLGAMCGAGL V 

TM4 



KAFQ- SAYF 
KAFQ- SAYF 
KGFQ- SAYY 
KGFQ- SAYY 
KGFQ- SAFY 
KGFQ- SGFY 
KGFQ- KSFY 
KAFHGAHWY 



DRYGGGANSL 
DRYGGGANSL 
VRYGGGANEL 
VRYGGGANEL 
VRYGGGANEL 
ARYGGGANEV 
NRYGGGVNTV 
LRYGGGANEL 

TM5 



ASGY 
ASGY 
SDGY 
SDGY 
SAGY 
SAGY 
SDGY 
AAGY 



SRGT 
SRGA 
SKGT 
SKGT 
SKGT 
STGT 
NKGT 
SKGA 



GLGAEI 
GLGAEI 
GLAAEI 
GL AAEI 
GLAAEI 
GLAAEI 
ALGAEI 
GLGAEI 



PVTGT 
PVTGT 
PI TGT 
PI TGT 
PI TGT 
PI TGT 
PVTGT 
PVTGT 



I GTFVLVYTVFSATDPK 
VGTFVL VYTVFSATDPK 
I GTFVLVYTVFSATDPK 
I GTFVLVYTVFSATDPK 
I GTFVLVYTVFSATDPK 
I GTFVLVYTVFSATDPK 
I GTFVLVYTVFSATDPK 
VGTFVL VYTVFSATDPK 

NPA LoopE 
G / NPARSLGAA y\ YNKD 
G l NPARSLGAA VVYNKD 
G l NPA ;SLGAA y\ YNKD 
YNKD 
YNND 



G l NPAHSLGAA \/\ 
G l NPARSLGAA \/\ 
G l NPARSLGAA VyYNHS 
G l NPARSFGPA \I\ FNND 
G l NPARSLGPA yyyHQR 



RNARDSHVPVLAPLPI GFAVFMVHLATI 
RNARDSHVPVLAPLPI GFAVFMVHLATI 
RSARDSHVPVLAPLPI GFAVFMVHLATI 
RSARDSHVPVLAPL PI GFAVFMVHLATI 
RNARDSHVPVLAPLPI GFAVFMVHLATI 
RNARDSHVPVLAPLPI GFAVFMVHLATI 
RNARDSHVPVLAPLPI GFAVFMVHLATI 
RKVRDSHVPVLAPLPI GFAVFMVHLATI 

TM6 

KPWDDHWI FWVGPLVGAAI AAFYHQYI LRAGAI KALGSFRSNA 
KPWDDHWI FWVGPLLGAAI AAFYHQYI LRAGAI KALGSFRSNA 
KAWDDQWl FWVGPLI GAAI A A A Y HQY VL R AS AT K- LGSYRSNA 
KAWDDQWI FWVGPLI GAAI AAAYHQYVL RAS AT K- LGSYRSNA 
KAWDDHWl FWVGPFI GAAI AAAYHQYVL RAS AAK- LGSSASFSR 
KAWSDQWI FWVGPFI GAAI AALYHQI VLRASARG- YGSFRSNA 
KAWDDQWl YWVGPFVGAAVAAI YHQYI LRGSAI KALGSFRSNA 
KAWEDHWI FWVGPLI GA AA A ML YHQL VL RAGAA K A F AS F RNNQH F 



98 
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98 

97 

93 

97 

95 

88 



148 
150 
148 
147 
143 
147 
145 
138 

197 
199 
197 
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247 

249 
247 
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290 
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289 
288 
285 
288 
287 
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Figure 2. Allgnrrient of ZmPIP2;8 and other ZmPIP2 sequences showing the conserved NPA domains of loop B In bold, loop E In bold italics 
and both underlined." The transmembrane domains are labeled TM1-6. The conserved H of loop D is starred. The conserved Ar/R filter is in red. 
ZmPIP2;1(AAK26758.1), ZmPIP2;2 (AAK26759.1), ZmPIP2;3 (AAK26760.1), ZmPIP2;4 (AAK26761.1), ZmPIP2;5 (AAD28761.1), ZmPIP2;6 (AAK26762.1), 
ZmPIP2;7 (AAK26763.1). 



A closer look at the conserved amino acids involved in sub- 
strate-specific transport by NIPs may suggest what this MIP may 
be transporting. NIPs can be subdivided into three groups based 
on protein sequence''"" and ZmNIPl;! is in group I along with 
the first described NIP, GmNod26 from soybean. First the four 
residues of the Ar/R filter are conserved between GmNod26 and 
ZmNIPl ;1 and second the Ar/R filter of NIPls are more simi- 
lar to GlpF (Table 2). In fact all the NIPls have the same Ar/R 
motifs except AtNIP3;l, which has an I at the second position, 
which is also an aliphatic amino acid (Table 2). The ZmNIPl;l 
protein sequence has a hybrid P1-P5 residues partly consisting of 



the water channeling aquaporin and the glycerol channeling aqua- 
porin (Table 2). All the NIPls have glycerol like residues at PI and 
P5, but the water channeling type at P2-P4. Froger et al."" men- 
tions this hybrid type in their analysis of NLMl, later renamed 
AtNIPl;l. All of the NIPls (Table 2) have invariant PI (F), P3 (A) 
and P4 (Y). At P2 all were S except OsNIPl;3, which was a T, still 
a neutral amino acid. P5 had the greatest amount of variability, V, 
I or L all aliphatic residues. AtNIPl;l, AtNIPl;2 and GmNOD26 
have been shown to transport glycerol (as reported in).'" However 
GmNOD26 also transports ammonia and formamide. AtNIPl;l 
also transports arsenite and AtNIPl;2 also transports H^O^.^" 
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ZmPIP2;3 
- ZmPIP2;4 

ZmPIP2;7 



■ ZmPIP2;5 

ZmPIP2;8 

ZmPIP2;1 

ZmPIP2;2 



ZmPIP2;6 

ZmPIPI ;3/1 ;4 



0.05 



Figure 3. ZmPIP2;8 clusters with other ZmPIP2s in a phylogenetic tree comparing ZmPIP2;1 (AAK26758.1), ZmPIP2;2 (AAK26759.1), ZmPIP2;3 
(AAK26760.1), ZmPIP2;4 (AAK26761.1), ZmPIP2;5 (AAD28761.1), ZmPIP2;6 (AAK26762.1), ZmPIP2;7 (AAK26763.1), ZmPIPl;3/l;4 (Q9AQU5.1). The maximum 
likelihood method was used as described in the materials and methods. 



Is there a physiological role for an increase ZmNIPl;! tran- 
script and the transport of any of these non-water compounds 
during herbivory? Of the neutral molecules listed above 1^,0^ has 
been linked to herbivory. It is induced systemically by wound- 
ing and acts as a second messenger for the systemic induction of 
defense genes in tomato.^' It has been shown to play a role in early 
response to wounding and other stresses in ArabidopsisP In fact 
several studies show that many genes induced by infestation with 
either phloem feeding or chewing insects are also induced by oxi- 
dative stress."^ Finally, the maize inbred line MP708, which is 
resistant to several lepidopteran larvae have higher levels of JA, 
as well as higher levels of several JA responsive defense genes, 
higher levels of transcript for the enzyme, NADPH oxidase that 
synthesizes H^O, and also higher levels of H^O^ upon herbiv- 
ory in comparison to a susceptible inbred line.^* Perhaps upon 
SCR infestation of maize ZmNIPl;! is induced to transport an 
increased amount of H^O^ within the plant. Further work will be 
needed to confirm this prediction. 

Hove and Bhave"* identified several signatures for Hj^i trans- 
porters within the Ar/R filter, P1-P5 residues and additional spec- 
ificity determining residues (SDP). The consensus signatures for 
Hp^ transporters are as follows: Ar/R filter (H/F/W), (I/H/V), 
(A/T/G), (V/R); P1-P5 residues (T/Q/F), (A/S), A, (Y/F), (W/I); 
SDPs {SI A), (A/G), (L/V), (A/V/F/L), (I/V), (H/I/Q/L), (A/V), 
P."* The Ar/R filter of the H^, transporters could contain any 
of the amino acid residues found in the NIPls, ZmPIP2;8 and 
ZmTIP2;2. However the P1-P5 residues found in the H^O^ 
transporters differs in one position in 6 of 10 NIPls and one in 
the ZmPIP2;8. The P5 consensus residue in H^O^ transporters is 
either W/I, which is found in only five of ten NIPls. Comparing 
the NIPls, ZmPIP2;8 and ZmTIP2;2 for the H^O, transporter 
SDPs shows the following. The consensus residues in the first, 
second, third, fifth and eighth positions are present in all the 
NIPl proteins examined as well as in ZmPIP2;8 and ZmTIP2;2 
(Table 2). The remaining three positions vary from the consensus 



signature for a H^O^ transporting MIP at only one residue of the 
SDPs in seven of the ten NIPls examined (Table 2). ZmPIP2;8, 
AtNIPl;l, AtNIPl;2 and AtNIP3;l contain all H,0, transport- 
ers consensus SDPs. While ZmTIP2;2 differs with the consensus 
SDPs at the fourth and sixth position. 

Since AtNIPl;2 has been shown to transport H^O^^" are there 
other NIPls with the conserved amino acids that predict they 
also transport H^O^? Comparison of the NIPls and ZmTIP2;2 
and ZmPIP2;8 (Table 2) shows that AtNIP3;l does not vary 
in any of these positions from any of the known H^O^ MIPs. 
AtNIPl;!, OsNIPl;!, OsNIPl;4 and ZmPIP2;8 varies in only 
one of these conserved amino acids. ZmNIPl;!, OsNIPl;2, 
OsNIPl;3, GmNOD26, AtNIP2;l and ZmTIP2;2 vary in only 
two of the amino acids. Perhaps this natural variation in MIPs 
could be used to test the in silico results identifying important 
conserved amino acids. 

Finally, this study also describes a new aquaporin that had 
been overlooked in previous maize aquaporin research, but is not 
unique to the maize variety used in this study, since it is present in 
the sequenced and annotated maize genome from the well stud- 
ied maize variety B73. It contains the typical structural aspects of 
an aquaporin with six trans-membrane domains, the conserved 
NPA motifs and Ar/R filter that have been shown to be critical 
for the water channel in other well characterized aquaporins (Fig. 
2, Table 1). The only difference in the conserved E loop between 
ZmPIP2s is an A vs. a P respectively in the second to last con- 
served portion of this loop of ZmPIP2;7 and ZmPIP2;8. Indeed 
ZmPIP2;8 confers robust water transport capacity in Xenopus 
oocyte swelling assay (Fig. 4). In conclusion we have, perhaps by 
serendipity, identified a new maize aquaporin ZmPIP2;8. 

Materials and Methods 

Aquaporin cloning and qRT-PCR. For root RNA isolation, 1 g 
of ground root tissue was added to 10 ml Trizol (Invitrogen) and 
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Figure 4. The water permeability of the oocytes microinjected with the cRNAs of ZmPIP2;8 and ZmTIP2;2, the cRNA of /.occor/afa/co/or aquaporin 
JQ585595 as the positive control and nuclease-free water as the negative control. Swelling assay was conducted in 0.2x hypotonic MBM after 48 h 
incubation in MBM solution at 18°C. Bar represents the mean of 10 replicates and error bar represents standard deviation. Letters indicate significant 
difference at p < 0.05. 



incubated at room temperature for 5 min with frequent vortex- 
ing. Then 1 ml of chloroform was added vortexed for 15 sec, 
incubated at room temperature for 1 min and vortexed for 15 sec. 
The phases were separated by centrifugation at 20,190 x g for 
10 min in a centrifuge. The top phase was collected and precipi- 
tated with 1 volume of isopropanol overnight at -80°C. The pre- 
cipitate was pelleted by centrifugation at 11,950 x g for 10 min. 
The pellet was air-dried and resuspended in 200 jjlI sterile water. 
The samples were then further extracted using RNeasy mini spin 
columns with on-column DNase treatment as described above. 
RNA concentration and relative purity was assessed by absor- 
bance at 260 and 280 nm with a spectrophotometer and integrity 
examined on a 1% (w/v) agarose gel with TBE buffer. 

Primers for cloning were selected from the following sequences 
available at MaizeGDB (www.maizegdb.org) for NIP1;1- 
GRMZM2G04l980_TO2, for TIP2;2 GRMZM2G056908_ 
TOl and PIP2-8 GRMZM2G432926_T01. The following 
primers were used to clone the aquaporins. TIP2-2-f CTC AGA 
GCC TCA GAG CGC GAG GCA AGT G, TIP2-2r-GGG AAA 
TGA ATA CAA GGA GGG AAT GGA AGA GAG G, PIP2 
f-GGA AGA GGG AGG GGG GAG GGT GTT, PIP2 r-GAG 
AGA GGA AAT GGA TGG AGG GGG AAA AT, NlPl-lf- 
GGA GGA GGG GAG GAA GGA GGA GTT GGG, NlPl-lr- 
GGG GGG GAG GAG GTA GAG GGA GAT G. The Titanium 
One-step RT-PGR kit (Glontech) was used to produce cDNA, 
which was ligated to pGEM-T according to the manufacturer's 
protocol (Promega). cDNAs were sequenced and used to design 
custom Taqman primers and probes (see Appendix S4; Applied 
Biosystems). For measuring transcript levels during infestation. 



cDNA was synthesized from maize root RNA and qRT-PGR was 
performed as described in Lawrence et al.^ 

Protein sequence analysis. The SGR infestation-repressed 
ZmPIP2;8 described in Lawrence et al.^ is annotated in the maize 
genome and the transcript can be found at MaizeGDB (www. 
maizegdb.org) as GRMZM2G432926_T01, which resides on 
chromosome 5 of the B73 RefGen version 2 genome between 
38,210,596-38,202,808. A scan of the available maize aquaporin 
literature did not recover this sequence. The novel maize PIP2 
contains an open reading frame, encoding a putative 283 amino 
acid protein. Since this is a newly described ZmPIP2, it was com- 
pared with other ZmPIP2 protein sequences, showing most simi- 
larity to ZmPIP2;2 at 86%. To determine whether it contained 
the conserved sequences novel to PIPs, two trans-membrane pro- 
tein prediction programs were used, TMHMM Version 2.0 at 
www.cbs.dtu.dk/services/TMHMM/^' and HMMTOP at www. 
enzim.hu/hmmtop.^' The NIPl proteins were aligned using 
GLUSTALW2 (www.ebi.ac.uk/Tools/msa/clustalw2). The posi- 
tion of the conserved amino acid motifs, Ar/R, P1-P5 and SDPs 
were determined by comparison to alignments of known H^O^ 
transporters from Hove and Bhave,"* which include AtNIPl;2. 
These motifs were ascertained for ZmTIP2;2 and ZmPIP2;8 by 
aligning to ZmNIPl;! and other ZmPIP2s. 

Phylogenetic analysis of ZmPIP2;8. Alignment of 
sequences to ZmPIP2 proteins was performed using Megalign 
from lasergene 8.0 (DNASTAR). The phylogenetic trees were 
inferred using the Maximum Likelihood method based on 
the JTT matrix-based model. Initial tree(s) for the heuristic 
search were obtained automatically by applying Neighbor-Join 
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Table 2. Comparison of conserved amino acids involved in pore formation of the ZmPIP2;8, ZmTIP2;2 and the NIPIs from rice, ma\ze, Arabidopsis and 
soybean to the prototypical water channeling protein AQP1 and the prototypical glycerol channeling protein GIpF 

Ar/R filter P1-P5 residues SDPs H202 



Protein ID 


TIV12 


TM5 


LE1 


LE2 


PI 


P2 


P3 


P4 


P5 


TM3 


TM3 


TM3 


TM3 


TMS 


HE 


TM6 


TM6 


ZmNIP1;1 


W 


V 


A 


R 


F 


S 


A 


Y 


V 


A 


A 


L 


M 


V 


! 


V 


P 


OsNIP1;1 


W 


V 


A 


R 


F 


S 


A 


Y 


1 


A 


A 


L 


M 


V 


1 


V 


P 


OsNIP1;2 


w 


V 


A 


R 


F 


s 


A 


Y 


V 


S 


A 


L 


V 


V 


V 


A 


P 1 


OsNIP1;3 


w 


V 


A 


R 


F 


T 


A 


Y 


1 


S 


A 


L 


M 


V 


1 


A 


P 


1 OsNIP1;4 


w 


V 


A 


R 


F 


s 


A 


Y 


1 


A 


A 


L 


L 


V 


L 


L 


P 


AtNIP1;1 


w 


V 


A 


R 


F 


s 


A 


Y 


L 


S 


A 


L 


L 


1 


1 


V 


P 


AtNIP1;2 


w 


V 


A 


R 


F 


s 


A 


Y 


1 


S 


A 


L 


L 


V 


L 


V 


P 


GmNOD26 


w 


V 


A 


R 


F 


s 


A 


Y 


L 


S 


A 


L 


L 




1 


L 


P 


AtNIP2;1 


w 


V 


A 


R 


F 


s 


A 


Y 


L 


S 


A 


L 


L 




1 


L 




AtNIP3;1 


w 


1 


A 


R 


F 


s 


A 


Y 


1 


S 


A 


L 


V 




L 


V 


P 


ZmPIP2;8 


F 


H 


T 


R 


H 


s 


A 


F 


W 


A 


G 


V 


F 




H 


V 


P 


ZmTIP2;2 


H 


1 


G 


R 


I 


s 


A 


Y 


W 


A 


A 


L 


y 




N 


V 


P 


AQP1 


F 


H 


C 


R 


T 


s 


A 


F 


W 


















GIpF 


W 


G 


F 


R 


Y 


D 


K 


P 


L 



















Conserved amino acids with AQP1 protein from humans (blue) and GIpF a glycerol channel protein from £ coli (red) are shown. The invariant R residue 
of LE2 is shown in black. Consensus signatures for H^O^ transporting aquaporins are as follows: Ar/R filter (H/F/W), (l/H/V), (A/T/G), (V/R); P1-P5 residues 
(T/Q/F), (A/S), A, (Y/F), (W/l); SDPs (S/A), (A/G), (L/V), (A/V/F/L), (l/V), (H/l/Q/L), (A/V), P.'" The 8-SDP consensus residues from the H^O^ transporters con- 
served in the NIPIs, ZmTIP2;2 and ZmPIP2;8 are shown (green). Residues not found in the signatures of H^Oj transporters are underlined. AtNIP1;2 has 
been shown to be a functional H^O^ transporter.^" 



and BioNJ algorithms to a matrix of pairwise distances esti- 
mated using a JTT model and then selecting the topology with 
superior log likelihood value. The tree is drawn to scale, with 
branch lengths measured in the number of substitutions per 
site. The analysis involved nine amino acid sequences. All posi- 
tions containing gaps and missing data were eliminated. There 
were a total of 270 positions in the final data set. Evolutionary 
analyses were conducted in MEGAS."** 

Subcloning ZmTIP2;2 and ZmPIP2;8 and testing in 
Xenopus oocyte swelling assay. The full-length cDNA of gene 
ZmPIP 2;8 ofZea mays was sub-cloned from pGEM vector into 
the restriction site of Spel of the expression vector pXT7 con- 
taining the 77 promoter and the 5' and 3' UTR of Xenopus lae- 
vis |3-globin gene (provided by Dr. Warren Gallin, University 
of Alberta). The full-length cDNA of gene ZmTIP2;2 of maize 
was sub-cloned into pXT7 at the restriction site of Spel and 
EcoRI. The expression vector containing the target gene in the 
77 orientation was determined by sequencing and was linear- 
ized at the site of Ndel located downstream of 3' UTR of Xenopus 
(3-globin gene. The linearized vector was used as the template for 
in vitro synthesis of capped RNA by using T7 RNA polymerase 
(mMESSAGE mMACHINE T7 kit, Ambion). The cRNA of 
aquaporin JQ585595 in Laccaria bicolorsuzxn UAMH8232 was 
used as the positive control and nuclease-free water as the nega- 
tive control. In Xenopus oocyte swelling assay, 10 ng of cRNA or 
nuclease-free water was microinjected into each healthy oocyte 
at Stage V-VI after being treated with 2 mg/ml collagenase 



in modified Earth's media for 2 h ([MBM] 88 mM NaCl, 1 
mM KCl, 2.4 mM NaHCOj, 10 mM Hepes-NaOH, 0.33 mM 
Ca[N03]„ 0.41 mM CaCl,, 0.82 mM MgSO^, Gentamicin 
sulfate 0.05 g/1. Penicillin G 0.1 g/1 Na pyruvate 2.5 mM, pH 
7.5) followed by a hypertonic solution of 100 mM potassium 
phosphate, pH 6.5 as described previously,^'''" using an auto- 
matic nanoliter injector (Nanoject II, Drummond Scientific). 
After incubation in 200 mosmol MBM in scintillation vials at 
18°C for 48 h, the injected oocyte was transferred into MBM 
in one well of a four-well Petri dish to take the initial image 
by an Olympus microscopy camera connected to an Olympus 
compound light microscope under 4x objective lens using the 
image software QCapture Pro 6.0 (Qlmaging) and then was 
transferred into D = 0.2 hypotonic MBM (40 mosmol) in the 
adjacent well and immediately the serial images were captured 
every 10 sec interval for 3 min, in order to track the changes in 
oocyte volume due to water influx. The diameter and surface 
area of oocytes were analyzed using Image] (V.1.44o, http:// 
imagej.nih.gov/ij/index.html). The initial transmembrane 
volume flux and osmotic water permeability co efficiency {P^ 
were calculated based on Zhang and Verkman^' to represent the 
water permeability of the oocytes injected with cRNAs of each 
aquaporin gene. Statistic analysis was conducted using Origin 
8.0 (OriginLab). 

Disclosure of Potential Conflicts of Interest 
No potential conflicts of interest were disclosed. 



www.landesbioscience.com 



Plant Signaling & Behavior 



624937-7 



References 

1. Hui D, Iqbal J, Lehmann K, Gase K, Saluz HP, 
Baldwin IT. Molecular interactions between the special- 
ist herbivore Manduca sexta (lepidoptera, sphingidae) 
and its natural host Nicoriana arrenuata: V. microar- 
ray analysis and further characterization of large-scale 
changes in herbivore-induced mRNAs. Plant Physiol 
2003; 131:1877-93; PMID: 12692347; http://dx.doi. 
org/10.1104/pp.l02.018176 

2. Reymond P, Bodenhausen N, Van Poecke RMP, 
Krishnamurthy V, Dicke M, Farmer EE. A conserved 
transcript pattern in response to a specialist and 
a generalist herbivore. Plant Cell 2004; 16:3132- 
47; PMID:15494554; http://dx.doi.org/10.1105/ 
tpc.104.026120 

3. De Vos M, Van Oosten VR, Van Poecke RMP, Van Pelt 
JA, Pozo MJ, Mueller MJ, et al. Signal signature and 
transcriptome changes of Arabidopsis during patho- 
gen and insect attack. Mol Plant Microbe Interact 
2005; 18:923-37; PMID:16167763; http://dx.doi. 
org/10.1094/MPMI-18-0923 

4. Ralph SG, Yueh H, Friedmann M, Aeschliman D, 
Zeznik JA, Nelson CC, et al. Conifer defence against 
insects: microarray gene expression profiling of Sitka 
spruce {Picea sitchensis) induced by mechanical wound- 
ing or feeding by spruce budworms {Choristoneura occi- 
dentalis) or white pine weevils {Pissodes strobi) reveals 
large-scale changes of the host transcriptome. Plant Cell 
Environ 2006; 29:1545-70; PMID:16898017; http:// 
dx.doi.org/10.1111/j.l365-3040.2006.01532.x 

5- Ehlting J, Chowrira SG, Mattheus N, Aeschliman DS, 
Arimura GI, Bohlmann J. Comparative transcriptome 
analysis of Arabidopsis thaliana infested by diamond 
back moth {Plutella xylostella) larvae reveals signatures 
of stress response, secondary metabolism, and signal- 
ling. BMC Genomics 2008; 9:154; PMID:18400103; 
http://dx.doi.org/10.1186/l471-2164-9-154 

6. Lawrence SD, Novak NG, Ju CJT, Cooke JEK. 
Examining the molecular interaction between pota- 
to {Solanum tuberosum) and Colorado potato beetle 
{Leptinotarsa decimlineatd) . Botany 2008; 86:1080-91; 
http://dx.doi.org/10.1139/B08-074 

7. Lawrence SD, Novak NG, Kayal WE, Ju CJT, Cooke 
JEK. Root herbivory: molecular analysis of the maize 
transcriptome upon infestation by Southern corn root- 
worm, Diabrotica undecimpunctata howardi. Physiol 
Plant 2012; 144:303-19; PMID:22172013; http:// 
dx.doi.org/lO.llll/j. 1399-3054.201 1.01557.x 

8. Forrest KL, Bhave M. Major intrinsic proteins (MIPs) 
in plants: a complex gene family with major impacts on 
plant phenotype. Funct Integr Genomics 2007; 7:263- 
89; PMID: 17562090; http://dx.doi.org/10.1007/ 
S10142-007-0049-4 

9. Maurel C, Verdoucq L, Luu DT, Santoni V Plant aqua- 
porins: membrane channels with multiple integrated 
functions. Annu Rev Plant Biol 2008; 59:595-624; 
PMID:18444909; http://dx.doi.org/10.ll46/annurev. 
arplant.59.032607.092734 



10. Schaeffer ML, Harper LC, Gardiner JM, Andorf CM, 
Campbell DA, Cannon EK, et al. MaizeGDB: cura- 
tion and outreach go hand-in-hand. Database (Oxford) 
2011; 2011:bar022; PMID:21624896; http://dx.doi. 
org/ 1 0. 1 093/database/bar022 

11. Chaumont F, Barrieu F, Wojcik E, Chrispeels MJ, Jung 
R. Aquaporins constitute a large and highly divergent 
protein family in maize. Plant Physiol 2001; 125:1206- 
15; PMID:1 1244102; http://dx.doi.org/10.1104/ 
pp.125.3.1206 

12. Bienert GP, BienertMD, JahnTP, Boutry M, Chaumont 
F. Solanaceae XIPs are plasma membrane aquaporins 
that facilitate the transport of many uncharged sub- 
strates. Plant J 2011; 66:306-17; PMID:2 124 1387; 
http://dx.doi.Org/10.llll/j.1365-313X.2011.04496.x 

13. Gupta AB, Sankararamakrishnan R. Genome-wide 
analysis of major intrinsic proteins in the tree plant 
Populus trichocarpa: characterization of XIP subfamily 
of aquaporins from evolutionary perspective. BMC 
Plant Biol 2009; 9:134; PMID:19930558; http:// 

dx.doi.org/10.1186/l471-2229-9-134 

14. Danielson JA, Johanson U. Unexpected complexity of 
the aquaporin gene family in the moss Physcomitrella 
patens. BMC Plant Biol 2008; 8:45; PMID:1 8430224; 
http://dx.doi.org/10.1186/l471-2229-8-45 

15. de Groot BL, Frigato T, Helms V, Grubmijller H. 
The mechanism of proton exclusion in the aqua- 
porin-1 water channel. J Mol Biol 2003; 333:279- 
93; PM1D:14529616; http://dx.doi.0rg/lO.lOl6/j. 
jmb.2003.08.003 

16. Froger A, Tallur B, Thomas D, Delamarche C. 
Prediction of functional residues in water channels 
and related proteins. Protein Sci 1 998; 7: 1458- 
68; PM1D:9655351; http://dx.doi.org/10.1002/ 
pro.5560070623 

17. Ankala A, Kelley RY, Rowe DE, Williams WP, Luthe 
DS. Foliar herbivory triggers local and long dis- 
tance defense responses in maize. Plant Sci 2013; 
199-200:103-12; PMID:23265323; http://dx.doi. 
org/10.10l6/j.plantsci.2012.09.017 

18. Hove RM, Bhave M. Plant aquaporins with non-aqua 
functions: deciphering the signature sequences. Plant 
Mol Biol 2011; 75:413-30; PMID:21308399; http:// 
dx.doi.org/10.1007/slll03-011-9737-5 

19. Liu Q, Wang H, Zhang Z, Wu J, Feng Y, Zhu Z. 
Divergence in function and expression of the NOD26- 
like intrinsic proteins in plants. BMC Genomics 
2009; 10:313; PMID:19604350; http://dx.doi. 
org/10.1186/1471-2164-10-313 

20. Dynowski M, Schaaf G, Loque D, Moran O, Ludewig 
U. Plant plasma membrane water channels conduct the 
signalling molecule Hp^. Biochem J 2008; 414:53- 
61; PMID:18462192; http://dx.doi.org/10.1042/ 
BJ20080287 

21. Orozco-Cardenas ML, Narvaez-Vasquez J, Ryan CA. 
Hydrogen peroxide acts as a second messenger for the 
induction of defense genes in tomato plants in response 
to wounding, systemin, and methyl jasmonate. Plant 
Cell 2001; 13:179-91; PMID:1 1 158538 



22. Miller G, Schlauch K, Tam R, Cortes D, Torres MA, 
Shulaev V, et al. The plant NADPH oxidase RBOHD 
mediates rapid systemic signaling in response to diverse 
stimuli. Sci Signal 2009; 2:ra45; PMID: 19690331; 
http://dx.doi.org/ 10.11 26/scisignal.2000448 

23. Kerchev PI, Fenton B, Foyer CH, Hancock RD. Plant 
responses to insect herbivory: interactions between 
photosynthesis, reactive oxygen species and hormonal 
signalling pathways. Plant Cell Environ 2012; 35:441- 
53; PMID:21752032; http://dx.doi.org/10. 1111/ 
j. 1365-3040.201 1.02399.x 

24. Shivaji R, Camas A, Ankala A, Engelberth J, Tumlinson 
JH, Williams WP, et al. Plants on constant alert: 
elevated levels of jasmonic acid and jasmonate-induced 
transcripts in caterpillar-resistant maize. J Chem Ecol 
2010; 36:179-91; PMID:20148356; http://dx.doi. 
org/1 0. 1 007/s 1 0886-0 1 0-9752-z 

25. Sonnhammer ELL, von Heijne G, Krogh A. A hidden 
Markov model for predicting transmembrane helices 
in protein sequences. In: Glasgow J, Littlejohn T, 
Major F, Lathrop R, Sankoff D, Sensen C, eds. Proc. of 
Sixth Int. Conf on Intelligent Systems for Molecular 
Biology Menio Park, CA: AAAI Press, 1998: 175-182. 

26. Tusnady GE, Simon 1. The HMMTOP transmem- 
brane topology prediction server. Bioinformatics 
2001; 17:849-50; PM1D:1 1590105; http://dx.doi. 
org/10.1093/bioinformatics/17.9.849 

27. Jones DT, Taylor WR, Thornton JM. The rapid genera- 
tion of mutation data matrices from protein sequences. 
Comput Appl Biosci 1992; 8:275-82; PMID:1633570 

28. Tamura K, Peterson D, Peterson N, Stecher G, Nei 
M, Kumar S. MEGA5: molecular evolutionary genet- 
ics analysis using maximum likelihood, evolutionary 
distance, and maximum parsimony methods. Mol 
Biol Evol 2011; 28:2731-9; PMID:21546353; http:// 
dx.doi.org/ 1 0. 1 093/molbev/msr 12 1 

29. Cao Y, Anderova M, Crawford NM, Schroeder JI. 
Expression of an outward-rectifying potassium chan- 
nel from maize mRNA and complementary RNA 
in Xenopus oocytes. Plant Cell 1992; 4:961-9; 
PMID:1392603 

30. Sand RM, Atherton DM, Spencer AN, Gallin WJ. 
jShawl, a low-threshold, fast-activating K(v)3 from the 
hydrozoan jellyfish Polyorchis penicillatus. J Exp Biol 
2011; 214:3124-37; PMID:21865525; http://dx.doi. 
org/10.1242/jeb.057000 

31. Zhang RB, Verkman AS. Water and urea perme- 
ability properties of Xenopus oocytes: expression of 
mRNA from toad urinary bladder. Am J Physiol 1991; 
260:C26-34; PMID:1987778 



e24937-8 



Plant Signaling & Behavior 



Volume 8 Issue 8 



